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Abstract 

The inner layers of a neutron star crust, composed of a Coulomb lattice of neutron rich 
nuclear clusters immersed in a sea of “free” superfluid neutrons, are closely analogous to 
periodic condensed matter systems such as electronic, photonic or phononic crystals. Applying 
methods from solid state physics to the neutron star context, we study the transport properties 
of those “free” neutrons for the outermost layers of the inner crust, near the drip point 
Pdrip ~ 4 X 10^^ g.cm“^. In particular, we evaluate the effective neutron mass resulting 
from Bragg scattering by a band structure calculation. Comparison is made with the case of 
electrons in solids. The observational consequences are briefly discussed. 


1 Introduction 

The outer layers of a neutron star are formed of a solid crust composed of a Coulomb lattice of 
very neutron rich nuclei immersed in a nearly uniform relativistic electron gas. In the deeper 
layers corresponding to densities beyond the drip threshold pdrip — 4 x 10^^ g.cm“^, nuclei 
are embedded in a sea of “free” neutrons which are expected to become superfluids in mature 
neutron stars whose temperature has dropped below the critical temperature for the onset 
of superfluidity (for a review of neutron star crust matter, see PP and references therein). 
Below the crust, at densities above ~ 10^^ g.cm“^ the nuclei merge into a uniform mixture of 
nucleons and leptons. 

The crust of a neutron star, which represents only about 1% of the mass of the star and 
10% of its radius, is nevertheless important to understand observational phenomena, such as 
for instance pulsar glitches. Pulsars are strongly magnetised rotating neutron stars whose 
period, ranging from milliseconds to a few seconds, very slowly increases due to the loss 
of energy by electromagnetic and gravitational radiation. However some pulsars have been 
observed to suddenly spin up. These so called “glitches” are characterised by a frequency 
jump, which varies from JH/H ~ 10“® (Crab) up to JH/H ~ 10“® (Vela), often accompanied 
by a permanent change in the slow down rate from | 10"®-10"'^ (Crab) to lO-^-lO"^ 
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(Vela). The relaxation following a glitch lasts from days to years. It was suggested by Baym 
et al. [2] soon after the discovery of the first pulsars that such very long timescales indicate 
the presence of superfluid components in the interior of the star. It is now widely accepted 
that pulsar glitches are due to a sudden tranfer of angular momentum from the faster rotating 
neutron superfluid to the solid crust. However the detailed mechanism of these glitches remain 
uncertain. In order to understand such phenomena, it is necessary to investigate the transport 
properties of the neutron superfluid in the inner crust |H1 li¬ 
lt has been recently shown that the effects of the nuclear clusters on the superfluid neutrons 
lead to a renormalisation of the neutron mass due to Bragg scattering [^, which gives rise 
to the entrainment effect according to which the momentum of the neutron superfluid is not 
aligned with the corresponding velocity |3]. This effective mass has been calculated in the 
very deep layers of the crust close to the interface with the liquid core for spherical nuclei 
in] and for the “pasta” phases where the crust is formed by a lattice of slab (“lasagna”) or 
cylinder (“spaguetti”) shaped nuclei |Sj. The effective mass has been found to be close to 
the bare neutron mass in the bottom of the crust while at lower densities the effective mass 
is strongly increased, reaching values as high as m* ~ 15mn at the baryon density 5 x 10^^ 
g.cm“^. The purpose of the present work is to extend this calculation to much lower densities 
near the neutron drip point. 

2 Band theory and effective mass 

Since the classical paper of Negele & Vautherin [7j (subsequently refered simply as N&V), 
quantum calculations in the inner crust have been carried out within the Wigner-Seitz (W- 
S) approximation according to which the crust is decomposed into a compact arrangement 
of independent spheres centered around each nuclear cluster. While this approach greatly 
simplifies the determination of the equilibrium structure of the crust and its static properties 
such as the equation of state for instance, it cannot be applied for the study of transport 
properties. Since the free neutrons in the inner crust of a neutron star are closely analogous 
to free electrons in metals, a natural step forward beyong the W-S approximation is to apply 
the well know band theory. 

This theory is not only restricted to the description of electrons in solids, but has been 
also applied to photonic and phononic crystals in which light (photons) and sound (phonons) 
respectively play the role of the electrons. Within the band theory, the crust is assumed to 
be a perfect lattice whose composition is known. The clusters forming the crust are treated 
classically (this approximation is justified since each cluster contains about one hundred nu¬ 
cleons or more) while the free neutrons are considered as independent quasiparticles in some 
effective mean field which takes into account the interactions with the other nucleons (in the 
superfluid phase or in the clusters). 

2.1 Effective neutron mass 

Subject to a more or less arbitrary convention specifying the density uj of “free” neutrons, 
the effective neutron mass can then be expressed as jS] 


Within the BCS approximation, the “mobility” coefficient KL is given by an integral over the 
first Brillouin zone (B-Z) [5] 
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where £k.a is the single particle energy of Bloch wave vector k and band index a, lS.p is the 
Fermi surface averaged value of the pairing gap and // is the Fermi energy. 

At the densities of the inner crust region near the drip point for which we are interested, 
the pairing gaps obtained in microscopic calculations of homogeneous pure neutron matter 
are very small compared to the typical Fermi energies jH] and since the mobility coefficient 
is determined mainly by the single particle energies £-ka^ we shall neglect the pairing for the 
calculation of the effective neutron mass. In this case the equation reduces to a Fermi surface 
integral |Hj 


/C = 
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The effective mass m*, as defined by o and m, has been refered in the solid state physics 
literature as an “optical” effective (electron) mass since it appears in the calculation of the 
real part of the dielectric constant of cubic metals cni, which can be written in the regime 
UJT l(r being the relaxation time) as 


£{U)} ~ 1 - + Winter 


(4) 


where Wp* is an effective plasma frequency defined in terms of the density n/ of “free” electrons 
and of their effective mass m* by 


Wp* = 




m* 


(5) 


and einter is the contribution due to interband transitions which is essentially constant for suf¬ 
ficiently low frequencies. It can be seen from o that this plasma frequency can be expressed 
as = 47re^/C and is therefore independent of the arbitrary convention used to specify which 
electrons are “free”. 

It should also be remarked that the electric conductivity a is related to the mobility 
coefficient JC introduced previously by the simple formula 


a = e^r/C, 


( 6 ) 


where e is the electron charge and r is the relaxation time. Unlike the effective electron 
mass m*, the electric conductivity is obviously independent of the convention adopted for 
defining which electrons are considered to be “free”. This point has been usually somehow 
obscured in the litterature by writing the conductivity in the form a = nfe^T/ni-i, of the 
classical Drude’s model. 

This effective mass can therefore be deduced experimentally by studying the optical prop¬ 
erties of solids or measuring the electric conductivity (or equivalently the resistivity p = l/cr) 
provided the relaxation time r is known. Values are typically in the range m* ~ 1 — 2me for 
most elements except for beryllium which has a very small effective mass mi,jme = 0.46 and 
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lead whose estimates of effective mass reach mi^lrrif. = 2.71 as reported in reference m (see 
appendix for some numerical values of rrii,). The density nj of “free” electrons is dehned by 
counting for each atom in the crystal the number of electrons in the outermost atomic shell. 
For instance monovalent metals (such as alkali or noble metals) have only one “conduction” 
or “free” electron per atom in the crystal therefore rif = 1/Vceii where Vceii is the volume of 
the W-S cell. 

The mobility coefficient m can also be written as the integral over the Fermi volume 



of the trace of the local effective mass tensor m defined by 

1 = (S) 

f^dkidkj' 

While the components of this tensor can be either positive or negative, it is clear from m 
and © that the “optical” effective mass dehned by m is always positive. 

Let us mention that the concept of the local effective mass tensor, which is well known in 
solid state physics, has been introduced in the context of neutron diffraction rather recently 
m- The effective mass for neutrons propagating along Bragg planes is found to be very small 
compared to the bare mass by several orders of magnitude as expected from the analogy with 
electrons in solids (the electron effective mass at the edge of the band gap in semiconductors 
is reduced by about two orders of magnitude m)- On the other hand, the “optical” effective 
mass m* as dehned by o is an average of the local effective mass © over all occupied states, 
not only those for which the Bloch wave vector lies on a Bragg plane, and can therefore be 
expected to be much larger (it has already been pointed out that the electron effective mass 
is indeed usually increased in metals, see also appendix). 

2.2 Band model of free neutrons 

The structure of the inner crust has been determined by N&V within a density matrix expan¬ 
sion method whose equations are similar to those obtained with effective two body nucleon- 
nucleon interactions of the Skyrme type in the Hartree-Fock approximation. We take their 
results for the lattice spacing and the composition of the nuclear clusters. The equilibrium 
structure of the crust is formed of a body centered cubic lattice. 

From the self-consistent calculations of N&V, we construct an effective single particle 
equation for the free neutrons, neglecting spin-orbit coupling as was done by N&V (this 
approximation can be justified from the fact that the nuclear clusters, which are very neutron 
rich, have a very diffuse surface) 


H = -V 


^ka^ka ) 


V + Un{r} . 


(9) 


( 10 ) 


2m^{r} 

As a result of the lattice symmetry, the single particle equations have to satisfy the Floquet- 
Bloch theorem m 

^ka{r + T} = (11) 
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where T is any lattice translation vector and k is a wave vector belonging to the first BZ. The 
equation can therefore be solved inside the W-S cell of the lattice (a truncated octahedron) 
with the boundary conditions CU). In the W-S approximation, this cell is replaced by a sphere 
with more or less arbitrary boundary conditions. For instance, N&V imposed the vanishing 
of the wave function or its radial derivative on the sphere depending on its parity. In this 
approximation, the k dependence of the single particle states (resulting from the translational 
symmetry) is therefore neglected. 

The kinetic term —7i^/2m®{r} and the mean field potential Uni^} are deduced from the 
self-consistent neutron and proton densities n„{r} and np{r} obtained by N&V, by applying 
the extended Thomas-Fermi approximation up to second order with Skyrme interactions m 


2m® {r} 


■^^ + ^(^i{2 + xi)+t2{2+x2)^nb{r}-^(^i{l + 2xi)-t2{l + 2x2)^nn{r} , ( 12 ) 


Un{r} = ^to [(2 -h xo)nb{r} - (I + 2xo)n„{r}] - [(2 -F X 3)(2 -h 'y)nb{ry^^ 

- (2X3 + l)(2nb{r}^nn{r} -h 'ynb{ry~^{np{r}‘^ + nn{r}^))] 

+ ^[ti (2 -f xi) +t2{2 + X2)]T{r} + ^[f2(2x2 + 1 ) - ti( 2 xi -b l)]r„{r} 

o o 

+ ^ [t2(2-b X 2 )-3fi(2-bxi)] V^nb{r} 

+ ^[3fi(2xi + l)+f2(2x2 + l)]vVW (13) 


where the kinetic densities are given by {q = n^p for neutron and proton respectively) 
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where is the well known Thomas-Fermi expression 
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and the second order contribution Tq = r®-b is given by the Bethe-Weiszacker term 
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plus non local terms arising from the Skyrme masses / = mjm® (calculated here for pure 
matter of the species q) 

2 


.r. _ 1 Vn, • V/ , 1 A/ 
/ e'"'' / 


6 


12 ’^'^ 


Y1 

f 


(17) 


We have used the SkM m parametrisation for the Skyrme force which is expected to give 
results close to those of N&V since its parameters like those in the density matrix expansion 
method of N&V, were adjusted on the neutron matter calculation of Siemens &: Pandharipande 
m and since this force was used a few years later by Bonche Vautherin m in the case of 
hot neutron star matter. For comparison, we have also used the more recent parametrisation 
SLy4 da Ca EHI, which has been specifically constructed for astrophysical purposes and 
which has been already applied in neutron star crust calculations mmiisiiiiiniin]- The 
parameters of these two Skyrme interactions are reported in the appendix. 
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3 Numerical methods 


3.1 LAPW method 


We have solved the equation Q by the standard solid state physics LAPW method . This 
method, whose idea was first introduced a long time ago by Slater m for the treatment 
of electrons in solids, provides a natural extension of the W-S approximation by combining 
the spherical symmetry aroung each nuclear cluster with the translational symmetry of the 
lattice. 

The W-S cell is decomposed into two regions: a sphere of radius Rg surrounding the cluster 
and an intersticial region, where in the muffin tin approximation the Skyrme mass and the 
mean field are assumed to take uniform values. 

The equation m is solved by expanding the wave function on a mixed basis set 


^k{r} = . 


(18) 


In the intersticial region the basis functions are defined by plane waves 


(pakir} = 


1 


J {k+Ka) r 


(19) 


where are reciprocal lattice vectors, which ensures that the Floquet-Bloch theorem m 
is satisfied. Inside the Slater sphere, the basis functions are given by 

^max ^ / \ 

(pakir] = EE AimUt{r] + Bimudf'} j yim{r} (20) 

i=0 Tn=—l ^ 


where Yim are spherical harmonics (r refers to the angular coordinates of r) and U£ is the radial 
solution of the equation inside the sphere for a given energy Ei- The dot means derivative 
with respect to the energy. The coefficients and are determined by imposing the 
continuity of the wave function and of its gradient on the sphere (see m for further details). 

From the Rayleigh-Ritz variational theorem, it follows that the Schrodinger equation © 
reduces to a generalised eigenvalue problem 


P P 

where Hap is the matrix of the Hamiltonian UlOj) 

Hap = {(pak\H(ppk ), 


( 21 ) 


( 22 ) 


and Sap is the overlap matrix 

Sap = {(pak\(ppk) ■ (23) 

The expansion m is truncated such as to include all reciprocal lattice vectors satisfying 


2mr 


^ -E^cutoff • 


(24) 


The choice of the radius Rg of the Slater sphere and of the linearisation energies {£i} is 
rather delicate due to the occurence of “ghost” bands. These bands can be however easily 
identified from their low k dependence (“flat” bands in the energy spectrum) and from their 
high sensitivity to small changes of the parameters Rg and (for a discussion on this issue, 
see for instance (23 )• 
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Figure 1; Fermi surface of copper obtained by Choy et al. jSl] (available online at 
http://www.phys. ufl. edu/fermisurface/). 

3.2 Fermi surface integration 

Fermi surface integrations were carried out within the Gilat-Raubenheimer scheme [^. The 
first B-Z is decomposed into small identical cubes, inside which the Fermi surface (if any) 
is approximated by a plane. From the knowledge of the energy bands and their gradient, it 
can be shown that the integral can be performed analytically inside each of those cubes. The 
integral over the whole Fermi surface is then obtained by summing the contributions of all 
the cubes. We have computed the energy gradient VjfcTfc by a finite difference scheme. It 
seems at first sight that the application of the Hellmann-Feynman theorem [.'12| . from which 
the gradient is obtainable by the formula 

Vfc<^fc = - TfcVfc'S'Q/g") , (25) 

a, 13 ^ ^ 

would lead to more accurate results. However this method, which we already applied in 
previous works Eiini, yield poor results in the present case because the errors in the wave 
functions obtained from the LAPW method are of order of {£ — £//, while the errors in the 
energy bands are only of order {£ — £// |HT| . 

We have tested the code by computing the energy bands and the effective mass m* in 
the empty lattice model m- For comparison, we have also computed the energy bands and 
“optical” effective electron mass of copper, which has been very often used as a benchmark 
for band structure calculation methods. 

3.3 Application to copper 

Copper is a noble metal whose electronic atomic configuration is [Hr]3d^^4s^ ([^^] denotes 
the electronic configuration of argon). A typical feature of such crystals is the mixing or 
“hybridization” of s and d atomic shells which leads to significant distortions of the Fermi 
surface from a sphere as shown on figure Q (in contrast to alkali metals for which the Fermi 
surface is nearly spherical). 

The closed atomic shells Is, 2s, 2p, 3s and 2>p give rise to fully occupied dispersionless 
energy bands in the crystal (“core” states) and it is therefore necessary to consider only the 
bands arising from the “valence” electrons in the 3d and 4s atomic shells. We have computed 
the energy bands of copper employing the classical Chodorov potential (see for instance issi 
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Figure 2: Band structure of copper (energies are shown in Rydbergs) along high symmetry 
lines of the first B-Z (with the conventional labelling HH])) computed with the LAPW method 
(see the text for further details). For comparison we have also plotted the energy spectrum 
(dashed line) of the ideal Fermi gas suitably translated into the first B-Z. The horizontal 
dotted line indicates the position of the Fermi energy. 


and references therein). Copper has a face centered cubic structure with a lattice spacing 
a = 3.6147 A. 

While the muffin-tin approximation is very good in neutron star crust, it is not so good 
in ordinary solids due to the long range of the Coulomb force. The radius Rg of the Slater 
sphere is thus set at the maximum possible value, namely half the nearest neighbor distance 
Rg = As already discussed by Koelling &: Arbman |31 j . the bands are rather sensitive 

to the choice of the linearisation energies Ei due to the presence of narrow (“localized”) d 
bands with low k dispersion. We thus calibrated the energies by computing the band 
structure with the original non linear APW method [28| . This method is more accurate 
but is computationally much more expensive than the linearised version. The expansion of 
the basis functions was truncated to Flcutoflf = 15 Ryds and I'max = 4. We found a fairly 
good agreement between the two methods with the following setting (energies are in Ryd): 
Sq = —0.90, Si = —0.20, and S 2 = S 3 = S^ = —0.70. As can be seen on figure |2l the energy 
bands show some resemblance to the ideal Fermi spectrum except for low k dispersion d bands 
in the range —0.55 < S < —0.8 Ryds. This is a characteristic feature of transition elements. 

We have used the Chodorov potential to compute the mobility coefficient of copper and the 
associated electron (“optical”) effective mass according to the equations o and ®. For this 
calculation we have fixed all the linearisation energies to the Fermi energy Si = = —0.385 

Rydbergs. In solid state physics, the number of “free” electrons is usually taken as the number 
of electrons in the “conduction” band. In the case of copper, there is one “conduction” electron 
per ion in the crystal from the atomic 4s shell and therefore nf = 1/Vceii- The corresponding 
(“optical”) effective electron mass as dehned by o and ©, m* ~ 1.285me obtained by 
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integration over 11375 cubes in the irreducible domain (therefore 546000 cubes in the B-Z), is 
in agreement with the earlier estimate of Seagall m for the same potential m*/me = 1.3±0.1. 
There is however some ambiguities in the specihcation of the “free” electron density nj m 
the case of transition metals since these elements may have different possible valences due to 
the presence of partially or completely filled d shell in the atom. We have also seen in the 
case of copper (see hgure|2) that in the crystal the electrons from the 3d atomic shell strongly 
hybridize with the “conduction” electrons from the outermost 4s shell. One may therefore 
adopt the convention of counting also these d electrons as “free”, in which case the effective 
mass becomes much larger m* ~ 14.14me since nf = 11/Vceii- 

4 Results and discussion 

We have considered the largest W-S cell of N&V (the radius of the W-S sphere is = 54 
fm) which corresponds to the onset of neutron drip at the baryon density p = 4.66 x 10^^ 
g.cm“^. The nucleon densities are expressed as Woods-Saxon functions inside the Slater 
sphere of radius Rs 


nq{r} = Uqo + {Uqi - Uqo)/ (^ + , (26) 

and are constants outside (and equal to nq{Rs}). The parameters were adjusted so as to fit 
the density profiles given by N&V, imposing that Upo = 0 and Uno = Pg where pc is the 
average density of the exterior neutron gas given by N&V. The lattice spacing a of the body 
centered cubic lattice was defined such that the volume of the W-S polyhedron, which is given 
by Vceii = a^/2 is equal to the volume of the W-S sphere, which leads to the formula 

a = 2(-K/3f/^R^, . (27) 

We found good convergence with the energy cutoff Tcutoff = 1 MeV and ^max = 6. We fixed 
the radius of the Slater sphere to Rs = 35 fm and all energy parameters for 0 < £ < £max 
to the mid value 0.15 MeV for the range 0 < T < 0.3 MeV (the energy origin is chosen as 
the constant value of the potential outside the muffin tin). We compared the band structure 
with that of the non linear APW method |28| and we found good agreement with this choice 
of parameters. 

The single particle energy spectrum for the “free” neutrons is shown on figure 01 for wave 
vectors k along specific symmetry lines of the first B-Z. The figure shows that unlike the 
band structure of copper computed in the previous section (see hgureO the neutron energy 
spectrum is actually very close to that of an ideal Fermi gas, except for avoided crossings 
around some specihc high symmetry points. The reason of this similarity lies in the so called 
Phillips-Kleinman cancellation theorem m- As a result of the Pauli exclusion principle, 
the orthogonalization of the “free” neutrons states to the “core” states leads to a repulsive 
non local and energy dependent potential which essentially vanishes outside nuclei and which 
partially compensates the strong attractive nuclear potential inside nuclei. In the present 
case, as the nuclei are very far apart {Rn,Rp 'C R^s) and since the nuclear interactions are 
short range the screening is almost complete. Consequently the effects of the nuclear clusters 
on the free neutron gas are negligible for bulk properties involving all free or valence states. 
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Figure 3: Energy bands (in MeV) of the “free” neutrons in neutron star crust (solid line) 
along high symmetry lines of the first B-Z (using the conventional labelling |29j i compared 
to the energy spectrum of the ideal Fermi gas (dashed line) suitably translated into the first 
B-Z for comparison. The band structure has been obtained with SkM Skyrme forces (details 
are explained in the text). 


This is illustrated on figure |1] for the total energy density e of the free neutrons given by 

e{nf} = r SM{S}dS , nj = T M{S}dS (28) 

Jo Jo 


where AA{T} is the density of single particle states defined by an integral over the constant 
energy surface 8 ^ = 8 


N{8} 



(29) 


It can be seen on figure |1] that the total energy density of the free neutrons does not depend 
neither on the chosen Skyrme interaction nor on the presence of the clusters and is actually 
very close to that of the ideal Fermi gas. We have computed the Fermi surface integral (OH 
for each energy on a grid of 20000 points in the range 0 < T < 0.3 MeV. Each of the Fermi 
surface integrations were performed with 9920 cubes in the irreducible domain, i.e. 476160 
cubes in the first B-Z. Increasing the number of cubes did not lead to any noticeable change 
on the figure. 

On the other hand, transport properties of the neutron gas depend only on the states 
lying on the Fermi surface whose shape is determined by the interactions between the clusters 
and the neutron gas (unlike its volume which is simply given by Vf = (27r)^nn) as in the 
ideal Fermi gas). As in the proceeding work |^, the density nj of free neutrons is defined 
by counting all single particle states whose energy is larger than the maximum value of the 
mean field potential As can be seen on figurelHlthe effective neutron mass m* deviates 
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Figure 4; Total energy density e (in MeV per fm^) of the free neutrons as a function of the 
free neutron density Uf (in fm“^) compared with that of the ideal Fermi gas. The arrow 
indicates the actual free neutron density in the W-S cell considered by Negele &: Vautherin. 


significantly from the “bare” neutron mass rUn (as in the case of electrons in copper previously 
discussed in section ESI), and exhibits large variations due to band effects, except at very low 
dripped neutron densities nj for which m* ~ m„. 

This can be understood from the consideration that the effects of the nuclear clusters are 
very small whenever the Fermi wave length is much larger than the lattice spacing, or in other 
words when the Fermi volume lies well inside the first B-Z. This means that the Fermi surface 
is nearly spherical for densities much lower than the density for which the Fermi sphere touches 
the B-Z zone boundaries. Such a situation occurs for alkali metals. For a body centered cubic 
lattice, this density threshold is given by n/ = -v/^vr/SVceii (at this density, the radius of the 
Fermi sphere equals the distance from the center of the first B-Z to the center of a face), 
which yields uj ~ 2.2 x 10“® fm“^ for = 54 fm. Distorsions of the Fermi sphere are 
expected to occur for slightly lower densities due to the formation of necks as in copper (see 
figure^. As shown on figure El above this density the Fermi surface area is reduced compared 
to the area of the corresponding sphere because by breaking the full translational symmetry 
(thereby raising some degeneracies), the presence of the nuclear lattice leads to the opening 
of holes on the Fermi surface (via Bragg scattering). In the limit for which the translational 
symmetry is completely lost (as in the case of isolated nuclei or in the W-S approximation), 
the Fermi surface disappears and its area vanishes which obviously implies that /C = 0. The 
mobility coefficient can also vanish whenever the energy spectrum exhibits a gap at the Fermi 
level. This situation is realised in ordinary insulators and in semiconductors for which the 
electric conductivity a = e^r/C = 0 since JC = 0 (at zero temperature). However this is still 
an open issue in the case of neutron star crust. 
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Figure 5: Effective neutron mass m*/mn as a function of the free neutron density nj (in 
fm“^). The arrow indicates the actual free neutron density in the W-S cell considered by 
Negele & Vautherin. 



Figure 6; Fermi surface area (in fm“^) as a function of the free neutron density nj (in fm“^) 
compared to the area of the corresponding Fermi sphere. The arrow indicates the actual free 
neutron density in the W-S cell considered by Negele & Vautherin. 
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5 Conclusion 


The band theory has been applied to investigate the transport properties of the neutron 
superfluid in the outermost layers of neutron star crust, at densities around 4 x 10^^ g.cm“^. 
We have described the free neutrons with an effective Schrodinger equation, involving a mean 
local mass m®{r} and a mean potential Unir} deduced from the results of Negele&Vautherin 
[7] in the extended Thomas-Fermi approximation with Skyrme nucleon-nucleon interactions. 
Unlike the calculations carried out by these authors, we have solved the equations beyond the 
Wigner-Seitz approximation by applying Bloch boundary conditions. 

It has been found that the effects of the nuclear clusters are very small when the neutron gas 
is very dilute, meaning that the Fermi wavelength is much larger than the lattice spacing. A 
similar case is observed in alkali metals. However at higher densities, the effective neutron mass 
m* (which is analog to the electron optical effective mass in solid state physics) is significantly 
increased compared to the “bare” mass due to Bragg scattering (as it is also observed for the 
case of electrons in ordinary metals) and fluctuates with the density as a result of band 
effects. The results also show significant differences with the choice of the parametrisation of 
the Skyrme interaction. A more accurate evaluation of the effective neutron mass in the crust 
therefore requires a fully self-consistent treatment in which both the equilibrium structure of 
the crust and the transport properties are computed within the same nuclear model. However 
the main conclusions drawn in this paper are not expected to be altered. 

The present calculations performed in the outermost layers of the inner crust combined 
with those carried out recently in the bottom layers [Hiini, suggest that the effective neutron 
mass m*/mn could take very large values in the middle layers of the crust where the size 
of the nuclear clusters is of the same order as the lattice spacing. This raises the question 
whether in some layers the effective neutron mass actually diverges (in the sense that the 
mobility vanishes /C = 0) due to the existence of band gaps in the energy spectrum. This 
question is one of the main issues in the study of any periodic materials such as for instance 
photonic or phononic crystals which have triggered considerable interests, both experimentally 
and theoretically. In our calculation, no such gaps were found for the outermost layers of the 
inner crust however as shown by Economou et al. isni in the context of photonic crystals, this 
does not exclude the possibility of gaps at higher densities (higher energies), in deeper layers. 
The existence of band gaps would have important implications for the dynamics of the star 
and for its cooling properties. 

Very large values of the effective mass have important consequences for the understanding 
of pulsar glitches. These sudden increases in the rotational frequency of some pulsars are 
interpreted as discontinuous transfers of angular momentum between the neutron superfluid 
and the solid crust. The present results suggest that in the middle layers the neutron superfluid 
would be strongly coupled to the crust so that such layers would contribute very little to the 
transfer of angular momentum. In the extreme case of band gaps, the superfluid would be 
locked to the solid parts so that relative currents could not develop (for the same reasons 
that electric currents do not flow in ordinary insulators). Large enough effective masses could 
also trigger a Kelvin-Helmholtz instability [H8| . which could provide a possible explanation 
for the origin of pulsar glitches. Further observational as well as theoretical investigations are 
required in order to sheld light on the nature of pulsar glitches and the link with the physics 
of “neutronic” crystals in neutron star crust. 
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Appendix 1: effective electron mass in metals 

Effective electron optical mass for a few elements taken from the compilation of reference El- 


metal 

m*/me 

Ag 

0.99-1.04 

A1 

1.52 

Au 

1.1 

Be 

0.46 

Ni 

1.7 

Ca 

1.8 

Cr 

1.72 

Fe 

1.9 

Cu 

1.45 

Ir 

1.41 

K 

1.08 

Li 

1.45 

Pb 

2.12-2.71 

Mo 

1.43-2 

Na 

0.98-1.25 

Nb 

1.86 

Pd 

1.66 

Rh 

1.87 

Ta 

1.65 

Ti 

1.39 

V 

1.63 

W 

1.31 


Appendix 2: LAPW matrix elements 

We have extended the LAPW method as discussed in the solid state context by Koelling and 
Arbman m in order to allow for a non uniform mass in the kinetic term. The end result is 
simply the replacement of the electron mass rrie by the Skyrme neutron mass in the 

expression of the matrix elements. In reference El, the authors used atomic units in which 
the kinetic factor /2me is equal to unity and therefore the electron mass does not appear 
explicity. We report in this appendix the full expressions for the matrix elements. 

The radial solution being normalised as 

rRs 

/ drr‘^Ui{r}‘^ = l, (30) 

Jo 
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the Hamiltonian and the overlap matrix elements are given by 




• q/3 (. 47ri?^ ji{|q« - q/3|} \ 

2mt{Rs}\°'^ Vcell |qa-q/3| ) 


( \ 

\2m®{Rs}) 


47r 

'hcell 


+ 00 


^{2£ + l)P£{cos e] £i{a(,{qo,]a(,{qp] + Neb(,{qf^]bi{qp]) + 7^^ 

(31) 


£=0 


{^(j)fy\(f)p^ — baf5 


A'KRlji{\cia - q/3|} 
Vcell |qa-q/3| 


+ 


47ri?^ 

Vcell 


+00 

E 


t=0 


{ai{qa]a^{q0] + N^bi{qa]hi,{q0]) 


(32) 


where is the regular spherical Bessel function of order i, Pf is the Legendre polynomial of 
order £, 6 is the angle between the vectors ^ k + and = k + K 0 (where and K 0 
are reciprocal lattice vectors), 



(33) 


aniq] = Ui{Rs}j[{qRs] - u'(,{Rs}3t{qRs] 


(34) 


h{q] = u'i{Rs}j^{qRs} - U(,{Rs}j[{qRs} 


(35) 


7^/3 = ui{Rs}u'(^{Rs}y'i{qaRs}jM0Rs} + j[{q0Rs}ji{qo,Rs]j 

- u'^{Rs}u'^{Rs}jl{qaRs}ji{q 0 Rs} - Ui{Rs}ui{Rs}j'ii{qaRs}j't{q 0 Rs} ■ (36) 
The prime means derivative with respect to the radial coordinate r. 


Appendix 3: Skyrme parametrizations 

The values of the parameters of the Skyrme forces used in this work are given in the following 
table (the unit of energy is MeV and the unit of length is fm) 



SkM 

SLy4 

^0 

-2645 

-2488.91 

tl 

385 

486.82 

t2 

-120 

-546.39 

h 

15595 

13777 

Xo 

0.09 

0.834 

Xl 

0 

-0.3438 

X2 

0 

-1 

X3 

0 

1.354 

7 

1/6 

1/6 
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